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ABSTRACT 

Using cosmological simulations, we address the interplay between structure and 
star formation in high-redshift galaxies via the evolution of surface density profiles. 
Our sample consists of 26 galaxies evolving in the redshift range z = 7—1, spanning the 
stellar mass range (0.2—6.4) x at z = 2. We recover the main trends by stacking 

the profiles in accordance to their evolution phases. Following a wet compaction event 
that typically occurs when the stellar mass is ~ 10®'® Mq at z ~ 2 — 4, the gas 
develops a cusp inside the effective radius, associated with a peak in star-formation 
rate (SFR). The SFR peak and the associated feedback, in the absence of further gas 
inflow to the centre, marks the onset of gas depletion from the central 1 kpc, leading to 
quenching of the central SFR. An extended, star-forming ring that forms by fresh gas 
during the central quenching process shows as a rising specific SFR (sSFR) profile, 
which is interpreted as inside-out quenching. Before quenching, the stellar density 
profile grows self-similarly, maintaining its log-log shape because the sSFR is similar 
at all radii. During the quenching process, the stellar density saturates to a constant 
value, especially in the inner 1 kpc. The stellar mass and SFR profiles deduced from 
observations show very similar shapes, consistent with the scenario of wet compaction 
leading to inside-out quenching and the subsequent saturation of a dense stellar core. 
We predict a cuspy gas profile during the blue nugget phase, and a gas-depleted core, 
sometimes surrounded by a ring, in the post-blue nugget phase. 

Key ■words: galaxies: evolution — galaxies: formation — galaxies: quenching — 
galaxies: high-redshift — galaxies: structure — galaxies: fundamental parameters 


1 INTRODUCTION 

At the present epoch, many early-type galaxies are massive 
systems, with stellar masses of M^, > 10^^ Mq, central stel¬ 
lar surface densities above 10^® Mq kpc“^, and hosting only 
little ongoing star-formation. In the high-redshift universe 
(z = 1 — 4), the progenitors of such massive galaxies had 
a significant gas-rich disc-like component and were forming 
stars at rates of up to a few hundred Mq yr“^. In this pa¬ 
per, we investigate the physical mechanism responsible for 
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the cessation of star formation, or “quenching”, and shed 
light on the interplay between the transformation of mor¬ 
phology and the changes in star-formation activity in these 
galaxies. 

There is solid observational evidence that the shut-down 
of star formation in galaxies is correlated with both galaxy 
mass and environment (e.g., Dressier 1980; Balogh et al. 
2004; Baldry et al. 2006; Kimm et al. 2009). The effects 
of galaxy stellar mass and environment in the quenching of 
galaxies appear to be separable up to redshift z = 1, i.e., it 
is possible to write the fraction of quiescent galaxies as the 
product of two functions, one a function of stellar mass only. 
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and the other a function of “environment” only (Peng et al. 
2010; Knobel et al. 2013; Kovac et al. 2014). Furthermore, 
quenching has been observed to correlate strongly with mor¬ 
phology and galaxy structure in the local universe (Kauff- 
mann et al. 2003; Franx et al. 2008; Robaina et al. 2012; 
Mendel et al. 2013; Cibinel et al. 2013a; Fang et al. 2013; 
Schawinski et al. 2014; Omand et al. 2014; Bluck et al. 2014; 
Woo et al. 2015) and at high -2 (Wuyts et al. 2011, 2012; 
Bell et al. 2012; Cheung et al. 2012; Szomoru et al. 2012; 
Barro et al. 2013; Lang et al. 2014; Tacchella et al. 2015c). 
On the other hand, quenched disc galaxies have also been 
found (McGrath et al. 2008; van Dokkum et al. 2008; Bundy 
et al. 2010; Salim et al. 2012; Bruce et al. 2012; Carollo 
et al. 2014). An important question is whether the separa¬ 
bility of the quenched fraction with mass and environment, 
highlighted by Peng et al. (2010), does imply two different 
quenching mechanisms, or is merely a reflection of the de¬ 
pendence on halo mass of the quenched fraction of satellite 
galaxies, as pointed out in Carollo et al. (2014). 

In this paper, we focus on quenching of galaxies at 
2 ~ 1 — 4. Even in the picture of two separate quenching 
mechanisms, at high redshifts the environmental quenching 
channel would be less relevant than the quenching channel 
that affects all galaxies independent of environment (Peng 
et al. 2010; Ilbert et al. 2013; Muzzin et al. 2013). Tacchella 
et al. (2015b, hereafter T15b) mapped out the M* and star- 
formation rate (SFR) distribution on scales of 1 kpc in 2 ~ 2 
star-forming galaxies (SFGs). They found that ~ 10^° Mq 
galaxies on the star-forming main sequence have fiat sSFR 
profiles, whereas ~ 10^^ Mq galaxies have sSFR profiles that 
rise toward the outskirts. Consistent results have been found 
at lower redshifts (e.g. Nelson et al. 2015). Furthermore, 
T15b used an empirical model to constrain the timescale 
quenching, finding that the massive 2 ~ 2 galaxies quench 
inside-out, where the star formation in the centre ceases 
within < 200 Myr, whereas the outskirts still form stars for 
1 — 3 Gyr. However, the physical nature of quenching, and its 
link to morphology, is yet to be understood. Several quench¬ 
ing mechanisms have been proposed. Galaxies can run out 
of gas by rapid gas consumption into stars, in combination 
with the associated outflows driven by stellar feedback (e.g., 
Dekel & Silk 1986; Murray et al. 2005) or super-massive 
black hole feedback (e.g., Di Matteo et al. 2005; Groton et al. 
2006; Ciotti & Ostriker 2007; Cattaneo et al. 2009), and/or 
by a slowdown of gas supply into the galaxy (Rees & Os¬ 
triker 1977; Dekel & Birnboim 2006; Feldmann et al. 2010; 
Hearin & Watson 2013; Feldmann & Mayer 2015). Morpho¬ 
logical / gravitational quenching suggests that the growth 
of a central mass concentration, i.e., a massive bulge, sta¬ 
bilizes a gas disc against fragmentation (Martig et al. 2009, 
2013; Genzel et al. 2014; Forbes et al. 2014). Long-term sup¬ 
pression of external gas supply happens once the halo mass 
grows above a threshold mass of ~ lO^^M©, either via virial 
shock heating (Birnboim & Dekel 2003; Keres et al. 2005; 
Dekel & Birnboim 2006; Keres et al. 2009), or by gravita¬ 
tional infall heating (Dekel & Birnboim 2008; Birnboim & 
Dekel 2011; Khochfar & Ostriker 2008), which can be aided 
by AGN feedback coupled to the hot halo gas (Dekel & Birn¬ 
boim 2006; Cattaneo et al. 2009; Fabian 2012). 

Regarding the morphology, the massive galaxies have 
transformed from a highly perturbed disc-like appearance 
in the rest-frame optical light at 2 = 1 — 4 to spheroidal-like 


appearance today. To understand how and when the galaxies 
assembled their stellar masses and developed structural com¬ 
ponents, the size-mass relation has been extensively studied 
at different redshifts (e.g., Daddi et al. 2005; Trujillo et al. 
2007; Cimatti et al. 2008; McGrath et al. 2008; van Dokkum 
et al. 2008; Szomoru et al. 2011; Newman et al. 2012; Barro 
et al. 2013; Dullo & Graham 2013; Poggianti et al. 2013; 
Shankar et al. 2013; Carollo et al. 2013b; Cassata et al. 
2013). Among other things, these studies found that the av¬ 
erage half-light radii of massive quiescent galaxies increase 
by a factor of 3 — 5 (at a constant mass, i.e., they were 
30 — 100-times denser) from redshift of 2 ~ 2 to the present. 
The origin of this evolution is debated. Often, this trend has 
been entirely ascribed to the physical growth of individual 
galaxies, where massive quiescent galaxies grow through ac¬ 
cretion of smaller galaxies in minor gas-poor mergers (e.g. 
Naab et al. 2009; Hopkins et al. 2009; Nipoti et al. 2009; 
Feldmann et al. 2010; Oser et al. 2012). On the other hand, 
several studies (Carollo et al. 2013b; Cassata et al. 2013; 
Poggianti et al. 2013) showed that the size growth of quies¬ 
cent galaxies can be explained by quenching of SFGs that 
keeps producing, at later epochs, new quenched galaxies 
with larger size than those of galaxies quenched at earlier 
epochs. The main challenge addressed here is to study the 
origin of a massive compact stellar core, with surface densi¬ 
ties reaching ~ 10^° M© kpc“^. Observations of 2 ~ 0 — 2.5 
galaxies inside radii of ~ 1 kpc indicate that these high cen¬ 
tral densities have formed early, before 2 ~ 2 (e.g., Wuyts 
et al. 2012; Saracco et al. 2012; Barro et al. 2013; Patel et al. 
2013; van Dokkum et al. 2010, 2013; Barro et al. 2014; van 
Dokkum et al. 2014; Morishita et al. 2015; Tacchella et al. 
2015b; van Dokkum et al. 2015). In this paper, we show that 
these densities must have originated from a dissipative com¬ 
paction of gas into a massive dense gas core that has largely 
converted into stars (Dekel & Burkert 2014). 

We know from cosmological simulations that high -2 
galaxies are continuously fed by streams from the cosmic 
web, consisting of smooth gas and a whole spectrum of 
merging clumps and galaxies (Birnboim & Dekel 2003; Keres 
et al. 2005; Dekel & Birnboim 2006; Ocvirk et al. 2008; Keres 
et al. 2009; Dekel et al. 2009a; Danovich et al. 2012, 2015). 
The high gas fraction and the high density of the universe at 
these high redshifts, combined with constant triggering by 
the intense instreaming of minor mergers, induce and main¬ 
tain violent disc instabilities (VDIs), which are characterized 
by turbulence and perturbations in the form of large tran¬ 
sient features and giant clumps (Noguchi 1998; Immeli et al. 
2004a,b; Genzel et al. 2006; Bournaud et al. 2007; Genzel 
et al. 2008; Agertz et al. 2009; Dekel et al. 2009b; Ceverino 
et al. 2010, 2012; Mandelker et al. 2014). To achieve the high 
gas mass in the centre, the gas has to flow to the centre at 
a rate that is faster than the SFR, such that the inflow is 
wet and the compaction is complete, not leaving behind an 
extended disc. 

Dekel & Burkert (2014) have introduced a toy model 
for describing the main aspects of this dissipative shrinkage 
of gaseous discs into compact star-forming systems. Zolotov 
et al. (2015, hereafter Z15) investigated the wet gas com¬ 
paction phase, and the way it leads to quenching, and iden¬ 
tified several physical mechanisms that can be responsible 
for the onset of compaction. They found that the wet com¬ 
paction is sometimes associated with VDI, and it is likely 
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triggered by an intense episode of gas inflow through gas- 
rich (mostly minor) mergers (e.g., Barnes & Hernquist 1991; 
Mihos & Hernquist 1996; Hopkins et al. 2006). The com¬ 
paction could alternatively be triggered by counter-rotating 
streams (Danovich et al. 2015), by low-angular-momentum 
recycled gas, and by tidal compression. Z15 put forward 
a picture in which galaxies evolve through compaction (to 
compact, star-forming galaxies, “blue nuggets”) and quench¬ 
ing phases. Following their definition of nuggets, we refer to 
a massive compact galaxy as a nugget, i.e., a galaxy that 
has a massive ~ 1 kpc core of high central density in stellar 
mass and possibly in gas density. Star-forming nuggets are 
called blue nuggets, but recall that the blue nuggets could 
actually be red due to dust. 

In a companion paper (Tacchella et al. 2015a, hereafter 
T15a), we have taken the analysis of Z15 one step further 
by investigating how the galaxies of these simulations evolve 
along the star-forming main sequence, and how they cease 
their star formation activity. We found that galaxies oscil¬ 
late about an equilibrium at the main-sequence ridge on 
timescales of ~ 0.4 tnubbie, ranging from 0.5 Gyr to 1 Gyr 
from 2 = 4 to 2. Galaxies on the upper part of the main 
sequence are compact, star-forming galaxies (blue nuggets), 
with higher gas fractions and shorter depletion times, con¬ 
sistent with observations (e.g., Magdis et al. 2012; Sargent 
et al. 2014; Huang & Kauffmann 2014; Genzel et al. 2015; 
Silverman et al. 2015; Scoville et al. 2015). We And that the 
blue nugget phase comes to an end when the balance of flow 
rates in the central 1 kpc tips to inflow rate < SFR -|- out¬ 
flow rate. The immediate reason is that while the SFR + 
outflow rate are at their maximum, the inflow is suppressed 
because the disc has largely disappeared (shrunk) and what¬ 
ever is left in the disc is low density. In addition, the bulge 
has grown, so the remaining disc is stable gravitationally 
(morphological quenching, see Martig et al. 2009; Genzel 
et al. 2014; Tacchella et al. 2015b) and does not drive more 
gas inward. If there is fresh gas supply from the halo on a 
sufficiently short timescale, it allows the disc to revive and 
a new compaction may be triggered to cause a new oscil¬ 
lation about the MS ridge. In order to maintain long-term 
quenching, one needs to prevent fresh gas supply from the 
halo. This is more likely at lower redshifts, where the cosmo¬ 
logical accretion timescale is longer than the depletion time, 
and especially so in a hot halo above a critical threshold halo 
mass (Birnboim & Dekel 2003; Keres et al. 2005; Dekel & 
Birnboim 2006). 

In this paper, we study the compaction and quench¬ 
ing of simulated galaxies using the time evolution of the 
density profiles of stellar mass, gas mass, and SFR. This is 
part of the effort to understand the interplay between the 
transformations of morphology and star-formation activity 
in massive, high -2 galaxies. We use high-resolution, zoom- 
in, hydro-cosmological simulations of 26 galaxies in the red- 
shift range 2 = 7 to 2 = 1. These are the same simulations 
as studied by Z15 and T15a. The suite of galaxies analyzed 
here were simulated at a maximum resolution of ~ 25 pc and 
include supernova and radiative stellar feedback. At 2 ~ 2, 
the halo masses are in the range Mvir ~ 10^^ — 10^^ Mq and 
the stellar masses are in the range (0.2 — 5.7) x 10^° M©. 

This paper is organized as follows. In Section 2 we give 
an overview of the simulations, including a summary of their 
limitations. In Section 3, we investigate surface density pro¬ 


files of individual galaxies and all galaxies based on average 
profles. Furthermore, we look into stellar mass and sSFR 
evolution in the centres as well as the outskirts of the simu¬ 
lated galaxies. Section 4 presents the emerging picture with 
the key figure showing average density profiles stacked by 
evolutionary stages of the galaxies. In Section 5 we present 
a comparison of the simulations with observational data, fo¬ 
cusing on the M* and SFR profile. Finally, we summarize 
our results in Section 6. 


2 SIMULATIONS 

To investigate the evolution of density profiles in high -2 
galaxies, we use zoom-in hydro-cosmological simulations of 
26 moderately massive galaxies, a subset of the 35-galaxy 
VELA simulation suite. Two papers (Ceverino et al. 2014; 
Z15) have presented the full details of these VELA simula¬ 
tions. Z15 and T15a used the same sample of 26 simulations 
and investigated similar questions concerning compaction 
and quenching. Z15 focused on the signatures of compaction 
in the simulation and T15a addressed the connection be¬ 
tween compaction, quenching and the confinement of the 
main sequence of SFGs. Additional analysis of the same suite 
of simulations are discussed in Moody et al. (2014), Snyder 
et al. (2015) and Ceverino et al. (2015b). In this section, we 
give an overview of the key aspects of the simulations and 
their limitations. 

2.1 Cosmological Simulations 

The VELA cosmological simulations make use of the Adap¬ 
tive Refinement Tree (ART) code (Kravtsov et al. 1997; 
Kravtsov 2003; Ceverino & Klypin 2009), which accurately 
follows the evolution of a gravitating N-body system and the 
Eulerian gas dynamics using an adaptive mesh refinement 
approach. The adaptive mesh refinement maximum resolu¬ 
tion is 17 — 35 pc at all times, which is achieved at densi¬ 
ties of ~ 10“"^ — 10® cm“®. In the circumgalactic medium 
(at the virial radius of the dark-matter halo), the median 
resolution amounts to ~ 500 pc. Beside gravity and hydro¬ 
dynamics, the code includes many physical process relevant 
for galaxy formation: gas cooling by atomic hydrogen and 
helium, metal and molecular hydrogen cooling, photoion¬ 
ization heating by the UV background with partial self¬ 
shielding, star formation, stellar mass loss, metal enrich¬ 
ment of the ISM and stellar feedback. Supernovae and stellar 
winds are implemented by local injection of thermal energy 
as in Ceverino & Klypin (2009); Ceverino et al. (2010) and 
Ceverino et al. (2012). Radiative stellar feedback is imple¬ 
mented at a moderate level, following Dekel & Krumholz 
(2013), as described in Ceverino et al. (2014). 

Cooling and heating rates are tabulated for a given gas 
density, temperature, metallicity and UV background based 
on the CLOUDY code (Ferland et al. 1998), assuming a 
slab of thickness 1 kpc. A uniform UV background based 
on the redshift-dependent Haardt & Madau (1996) model 
is assumed, except at gas densities higher than 0.1 cm“®, 
where a substantially suppressed UV background is used 
(5.9 X 10® erg s“® cm“® Hz“®) in order to mimic the partial 
self-shielding of dense gas, allowing dense gas to cool down 
to temperatures of ~ 300 K. The assumed equation of state 
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is that of an ideal mono-atomic gas. Artificial fragmentation 
on the cell size is prevented by introducing a pressure floor, 
which ensures that the Jeans scale is resolved by at least 7 
cells (see Ceverino et al. 2010). 

Star formation is assumed to occur at densities above a 
threshold of 1 cm“® and at temperatures below 10'* K. Most 
stars (> 90 %) form at temperatures well below 10^ K, and 
more than half of the stars form at 300 K in cells where the 
gas density is higher than 10 cm“^. The code implements a 
stochastic star-formation model that yields a star-formation 
efficiency per free-fall time of ~ 2 %, using the stellar ini¬ 
tial mass function of Chabrier (2003). At the given resolu¬ 
tion, this efficiency roughly mimics the empirical Kennicutt- 
Schmidt law (Kennicutt 1998). 

The thermal stellar feedback model incorporated in the 
code releases energy from stellar winds and supernova ex¬ 
plosions as a constant heating rate over 40 Myr following 
star formation. The heating rate due to feedback may or 
may not overcome the cooling rate, depending on the gas 
conditions in the star-forming regions (Dekel & Silk 1986; 
Ceverino & Klypin 2009). There is no artificial shutdown 
of cooling implemented in these simulations. Runaway stars 
are included by applying a velocity kick of ~ 10 km s“^ to 
30 % of the newly formed stellar particles. The code also 
includes the later effects of Type la supernova and stel¬ 
lar mass loss, and it follows the metal enrichment of the 
ISM. We incorporate radiation pressure through the addi¬ 
tion of a non-thermal pressure term to the total gas pres¬ 
sure in regions where ionizing photons from massive stars 
are produced and may be trapped. This ionizing radiation 
injects momentum in the cells neighbouring massive star 
particles younger than 5 Myr, and whose column density ex¬ 
ceeds 10^^ cm“^, isotropically pressurizing the star-forming 
regions. Ceverino et al. (2014) provides more details. 

The initial conditions for the simulations are based on 
dark-matter haloes that were drawn from dissipationless 
N-body simulations at lower resolution in three comoving 
cosmological boxes (box-sizes of 10, 20, and 40 Mpc/h). 
We assume the standard ACDM cosmological model with 
the WMAP5 values of the cosmological parameters, namely 
nm = 0.27, Da = 0.73, fib = 0.045, h = 0.7 and erg = 0.82 
(Komatsu et al. 2009). We selected each halo to have a given 
virial mass at 2 = 1 and no ongoing major merger at 2 = 1. 
This latter criterion eliminates less than 10 % of the haloes 
which tend to be in a dense, proto-cluster environment at 
2 ~ 1. The virial masses at 2 = 1 were chosen to be in the 
range Mvir = 2 x 10^^ — 2 x 10^^ Mq, about a median of 
4.6 X 10^^ Mq. If left in isolation, the median mass at 2 = 0 
was intended to be ~ 10^^ Mq. Realistically, the actual mass 
range is broader, with some of the haloes merging into more 
massive haloes that eventually host groups at 2 = 0. 

2.2 Limitations of the Current Simulations 

The cosmological simulations used in this paper are state-of- 
the-art in terms of high-resolution adaptive mesh refinement 
hydrodynamics and the treatment of key physical processes 
at the subgrid level. In particular, they trace the cosmolog¬ 
ical streams that feed galaxies at high redshift, including 
mergers and smooth flows, and they resolve the VDI that 
governs high -2 disc evolution and bulge formation (Ceverino 
et al. 2010, 2012, 2015a; Mandelker et al. 2014). 


Like in other simulations, the treatment of star forma¬ 
tion and feedback processes can still be improved. Currently, 
the code assumes a SFR efficiency per free fall time that is 
rather realistic, but it does not yet follow in detail the for¬ 
mation of molecules and the effect of metallicity on SFR 
(Krumholz & Dekel 2012). Furthermore, the resolution does 
not allow the capture of the Sedov-Taylor adiabatic phase 
of supernova feedback. The radiative stellar feedback as¬ 
sumed no infrared trapping, in the spirit of low trapping 
advocate by Dekel & Krumholz (2013) based on Krumholz 
& Thompson (2012). Other works assume more significant 
trapping (Murray et al. 2010; Krumholz & Dekel 2010; Hop¬ 
kins et al. 2012), which makes the assumed strength of the 
radiative stellar feedback lower than in other simulations. Fi¬ 
nally, AGN feedback, and feedback associated with cosmic 
rays and magnetic fields, are not yet implemented. Never¬ 
theless, as shown in Ceverino et al. (2014), the star forma¬ 
tion rates, gas fractions, and stellar-to-halo mass ratio are 
all in the ballpark of the estimates deduced from observa¬ 
tions, providing a better match to observations than earlier 
simulations. 

The uncertainties and any possible remaining mis¬ 
matches by a factor of order 2 are comparable to the ob¬ 
servational uncertainties. For example, the stellar-to-halo 
mass ratio is not well constrained observationally at 2 > 1. 
Recent estimates by Burkert et al. (2015) (see their Fig. 
5), based on the observed kinematics of 2 ~ 0.6 — 2.8 
SFGs, reveal significantly larger ratios than the estimates 
based on abundance matching (Conroy & Wechsler 2009; 
Moster et al. 2010, 2013; Behroozi et al. 2010, 2013a) at 
Mvir < Mq. We compare in Figure Al the stellar-to- 

halo mass ratio of our simulations with the ones from obser¬ 
vations. At 2 = 1, the simulated galaxies have a median halo 
mass of logjg Mvir = 11.7 ± 0.3 and a stellar to halo mass 
ratio of log^g M*/Mvir = —1.3 ± 0.2. Burkert et al. (2015) 
found logj^Q M»/Mvir = —1.5 ± 0.3 at logjg Mvir = 11.7, 
which is consistent with our simulations, while Behroozi 
et al. (2013a) found logj^g M*/Mvir = —1.8 ± 0.2 via abun¬ 
dance matching, which is a factor of 3 lower than our sim¬ 
ulated estimate. At 2 = 2, we find only little evolution in 
the halo-to-stellar mass ratio: our simulated galaxies have a 
median halo mass of log^g Mvu = 11-5 ± 0.3 and a stel¬ 
lar to halo mass ratio of logj^g M*/Mvir = —1.4 ± 0.2. 
Via abundance matching, Behroozi et al. (2013a) found 
logj^Q M*/Mvir = —2.3±0.2 at logjg Mvir = 11.5, but Burk¬ 
ert et al. (2015) found logj^g M*/Mvir = —1.6±0.3, which is 
again consistent with our simulations. Altogether, we con¬ 
clude that our simulations produce stellar-to-halo mass ra¬ 
tios that are in the ballpark of the values estimated from ob¬ 
servations, and within the observational uncertainties, which 
are also comparable to the uncertainties associated with the 
feedback recipes in the simulations. Therefore, it is sensible 
and adequate to use our simulations with the adopted feed¬ 
back prescription, while bearing in mind the factor-of-two 
uncertainties. 

It seems that in the current simulations, the compaction 
and the subsequent onset of quenching occur at cosmological 
times that are consistent with observations (see Fig. 12 of 
Zolotov et al. 2015 and Fig. 2 of Barro et al. 2013). However, 
with some of the feedback mechanisms not yet incorporated 
(e.g., fully resolved supernova feedback and AGN feedback), 
full quenching to very low sSFR values may not be entirely 
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reached in many galaxies by the end of the simulations at 
2 ; ~ 1. One should be aware of this limitation of the current 
simulations. Nevertheless, for all the above reasons, we adopt 
in this work the hypothesis that the simulations grasp the 
qualitative features of the main physical mechanisms that 
govern galaxy evolution through the processes of compaction 
and subsequent quenching. 

2.3 Simulated Galaxy Sample and Measurements 

From the suite of 35 galaxies, we have excluded three low- 
mass galaxies as well as six galaxies which have not been 
simulated down to 2 ; = 2.0. Therefore, our final sample con¬ 
sists of 26 galaxies. The virial and stellar properties of these 
galaxies are listed in Table 1. The virial mass Mvir is the to¬ 
tal mass within a sphere of radius i?vir that encompasses an 
overdensity of A( 2 :) = — 82 flA( 2 ) — 39 f 2 A( 2 :)^)/flm(z), 

where flA(z) and Q.rn.{z) are the cosmological parameters at 
2 ; (Bryan & Norman 1998; Dekel & Birnboim 2006). The 
stellar mass M* is measured within a radius of 0.2 x Rviv 

We start the analysis at the cosmological time corre¬ 
sponding to expansion factor a = 0.125 (redshift z = 7). At 
earlier times, the fixed resolution scale typically corresponds 
to a non-negligible fraction of the galaxy size. As can be seen 
in Table 1, most galaxies reach a = 0.50 {z = 1). Each galaxy 
is analysed at output times separated by a constant inter¬ 
val in a, Aa = 0.01, corresponding at 2 ; = 2 to ~ 100 Myr 
(roughly half an orbital time at the disc edge). Our sam¬ 
ple consists of totally ~ 900 snapshots in the redshift range 
2 = 6 — 1 from 26 galaxies that at 2 = 2 span the stellar 
mass range (0.2 —6.4) x 10^^ Mq. The half-mass sizes are 
determined from the M* that are measured within a radius 
of 0.2 i?vir and they range Ra — 0.4 — 3.2 kpc at 2 = 2. 

The SFR for a simulated galaxy is obtained by SFR = 

(Af*(tage tmax)/tniax)ta,ax 5 where A4* (tage ^ tmax) is the 

mass at birth in stars younger than tmax- The average 
(■)tmax is obtained by averaging over all fmax in the interval 
[40, 80] Myr in steps of 0.2 Myr. The fmax in this range are 
long enough to ensure good statistics. Our 26 galaxies have 
SFR ranging from ~ 1 to 33 Mq yr“^ at 2 ~ 2. 

The instantaneous mass of each star particle is derived 
from its initial mass at birth and its age using a fitting for¬ 
mula for the mass loss from the stellar population repre¬ 
sented by the star particle, according to which 10 %, 20 % 
and 30% of the mass is lost after 30 Myr, 260 Myr, and 2 
Gyr from birth, respectively. We consistently use here the 
instantaneous stellar mass, M*, and define the specific SFR 

by sSFR = SFR/MC. 

We split the simulated galaxies in a low- and high-mass 
sample based on the total stellar mass at 2 = 2. We divide 
the sample in about half at logj^Q M-^/Mq = 10.2. Galax¬ 
ies below and above this mass scale have a mean mass of 
logio M*/Mo = 9.9 ± 0.2 and log^o M^/Mq = 10.5 ± 0.2, 
respectively. This shows that our galaxies are relatively low- 


^ Note that this sSFR is not strictly the inverse of the timescale 
for the formation of the stellar system, for which the instanta¬ 
neous masses of stars have to be replaced by their initial masses. 
This makes only a small difference for our purpose here, given 
that in our galaxies the mass-loss rate is small compared to the 
SFR. 


and high-mass galaxies, and we do not study the most mas¬ 
sive or the least massive galaxies at 2 ~ 2 — 3. We caution 
that this is a rough division and that the galaxies which are 
in a certain mass bin at 2 = 2 do not necessarily belong to 
the same mass bin at all times. Generally, we find that the 
high-mass galaxies evolve through similar phases as the low- 
mass galaxies, but at earlier times and in a more pronounced 
way. 

The density profiles presented here are two dimensional 
surface density profiles E(r). We have chosen to present the 
surface density profiles instead of the three dimensional vol¬ 
ume density profiles p{r) because the surface density profiles 
are straight-forwardly accessible observationally. The deter¬ 
mination of the centre of the galaxy is outlined in detail 
in Appendix B of Mandelker et al. (2014). Briefly, starting 
form the most bound star, the centre is refined iteratively by 
calculating the centre of mass of stellar particles in spheres 
of decreasing radii, updating the centre and decreasing the 
radius at each iteration. We begin with an initial radius of 
600 pc, and decrease the radius by a factor of 1.1 at each 
iteration. The iteration terminates when the radius reaches 
130 pc or when the number of stellar particles in the sphere 
drops below 20. The three dimensional density profiles p(r) 
are then constructed from concentric spheres around this 
centre. We have used the simple relation E(r) = 2 • r • p(r) 
to convert the three dimensional density profiles to two di¬ 
mensional surface density profiles, which is a good approxi¬ 
mation for the face-on projection. The differences in surface 
density at a given radius for the face-on projection and our 
approach are < 20 %, i.e. well within our systematic uncer¬ 
tainty of the simulation. 


3 GROWTH OF GALAXIES 

We investigate how the surface density profiles of stellar 
mass, gas mass, SFR, and sSFR are evolving with time for 
individual simulated galaxies as well as for the whole sam¬ 
ple of galaxies. In particular, we use the evolution of the gas 
profile to study the wet compaction process and the sub¬ 
sequent phase of gas depletion from the core. We study the 
emergence of a dense stellar system through the evolution of 
the stellar density in the central region in comparison to the 
growth of the total mass. We then analyse the SFR density 
profile and the progression of quenching within the galaxies. 
Finally, we study the evolution of the half-mass sizes as a 
function of stellar mass and cosmic time. 


3.1 Evolution of the Surface Density Profiles 

Figure 1 shows the cosmic evolution of the stacked surface 
density profiles of stellar mass, gas mass, SFR, and sSFR for 
low-mass (left) and high-mass (right) galaxies, respectively. 
We split the simulated galaxies in a low- and high-mass sam¬ 
ple based on the total stellar mass at 2 = 2 , as described 
in Section 2.3. The individual profiles of the simulations are 
scaled before stacking to have the same effective density, 
Eefr = Qtot/{2'ivR^), where Qtot is the total relevant quan¬ 
tity and Re is the half-mass radius in 3D. The final, stacked 
profile has, by construction, a Eeg value that is equal to the 
median Eeff of all the galaxies in the sample. A complete 
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redshift 


Figure 1. Evolution of the stacked density profiles. From top to bottom, the surface density profiles of stellar mass (Em)? cold (T < 10"^ K) 
gas mass (Egas), SFR (Esfr)? s,nd sSFR are shown. The left and right panels show the stacks for the low- and high-mass galaxies, which 
are selected based on the stellar mass at z = 2. The colour coding corresponds to redshift, according to the colour bar at the bottom. 
The key feature is the development of a gas cusp in the massive galaxies that is associated with the wet compaction. Subsequently, the 
gas density in the centre decreases. The sSFR profile is fiat which leads to a self-similar growth of the stellar mass profiles. Towards 
lower redshifts, the sSFR decreases from inside out. The compaction with the succeeding central depletion occurs in different galaxies at 
different times, which smears out the evolutionary features in these stacks, especially for low-mass galaxies. The profiles for all galaxies 
are shown individually in Figures B1-B4. 
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The suite of 26 simulated galaxies. 


Galaxy 

Mvir 

IOI 2 Mq 
(2 = 2) 

M* 

10l° Mq 
(2 = 2) 

Mg 

10l° Mq 

(2 = 2 ) 

SFR 

Mq/jv 

(2 = 2 ) 

sSFR 
Gyr-i 
(2 = 2 ) 

-^vir 
kpc 
(^ = 2 ) 

Re 
kpc 
(2 = 2 ) 

^fin 

• 2 ^ fin 

VOl 

0.16 

0.22 

0.12 

2.65 

1.20 

58.25 

1.08 

0.50 

1.00 

V02 

0.13 

0.19 

0.16 

1.84 

0.94 

54.50 

2.23 

0.50 

1.00 

VOS 

0.14 

0.43 

0.10 

3.76 

0.87 

55.50 

1.73 

0.50 

1.00 

V06 

0.55 

2.22 

0.33 

20.72 

0.93 

88.25 

1.11 

0.37 

1.70 

V07 

0.90 

6.37 

1.42 

26.75 

0.42 

104.25 

3.41 

0.50 

1.00 

VOS 

0.28 

0.36 

0.19 

5.76 

1.58 

70.50 

0.81 

0.50 

1.00 

V09 

0.27 

1.07 

0.31 

3.97 

0.37 

70.50 

1.84 

0.39 

1.56 

VIO 

0.13 

0.64 

0.11 

3.27 

0.51 

55.25 

0.54 

0.50 

1.00 

Vll 

0.27 

1.02 

0.58 

17.33 

1.69 

69.50 

3.76 

0.46 

1.17 

V12 

0.27 

2.06 

0.19 

2.91 

0.14 

69.50 

1.27 

0.39 

1.56 

VIS 

0.31 

0.96 

0.98 

21.23 

2.21 

72.50 

5.36 

0.39 

1.56 

V14 

0.36 

1.40 

0.59 

27.61 

1.97 

76.50 

0.41 

0.42 

1.38 

V15 

0.12 

0.56 

0.14 

1.71 

0.30 

53.25 

1.33 

0.50 

1.00 

V20 

0.53 

3.92 

0.48 

7.27 

0.19 

87.50 

1.99 

0.44 

1.27 

V21 

0.62 

4.28 

0.57 

9.76 

0.23 

92.25 

1.89 

0.50 

1.00 

V22 

0.49 

4.57 

0.21 

12.05 

0.26 

85.50 

1.36 

0.50 

1.00 

V2S 

0.15 

0.84 

0.19 

3.32 

0.39 

57.00 

1.40 

0.50 

1.00 

V24 

0.28 

0.95 

0.28 

4.39 

0.46 

70.25 

1.87 

0.48 

1.08 

V25 

0.22 

0.76 

0.08 

2.35 

0.31 

65.00 

0.88 

0.50 

1.00 

V26 

0.36 

1.63 

0.25 

9.76 

0.60 

76.75 

0.80 

0.50 

1.00 

V27 

0.33 

0.90 

0.52 

8.75 

0.97 

75.50 

2.98 

0.50 

1.00 

V29 

0.52 

2.67 

0.39 

18.74 

0.70 

89.25 

2.47 

0.50 

1.00 

V30 

0.31 

1.71 

0.41 

3.84 

0.22 

73.25 

1.64 

0.34 

1.94 

V32 

0.59 

2.74 

0.37 

15.04 

0.55 

90.50 

2.66 

0.33 

2.03 

V33 

0.83 

5.17 

0.45 

33.01 

0.64 

101.25 

1.42 

0.39 

1.56 

V34 

0.52 

1.73 

0.42 

14.79 

0.85 

86.50 

2.10 

0.35 

1.86 


Table 1. Quoted are the total virial mass, Mvir, the stellar mass, M*, the gas mass. Mg, the star formation rate, SFR, the specific star 
formation rate, sSFR, the virial radius, Rvin the effective stellar (half-mass) radius, Rc (all at 2 = 2) and the final simulation snapshot, 
afin, and redshift, za^. The M*, Mg, SFR, and sSFR are measured within a radius of 0.2 X Rvir- 


list of all the surface density profiles for all galaxies individ¬ 
ually, ordered by total stellar mass at 2 = 2, can be found 
in Appendix B (Figures B1-B4). 

The evolution of the high-mass galaxies is quiet syn¬ 
chronized in cosmic time. Therefore, the key evolution¬ 
ary phases of compaction and quenching are visible in the 
stacked profiles (Figure 1, right panel). The gas density 
evolves from a flat profile at 2 > 4, through a cuspy pro- 
Hle at 2 ~ 3 — 4, to a ring-like gas distribution surrounding 
a gas-poor core at 2 < 2 — 3. These evolutionary stages are 
more noticeable in the individual profiles as these charac¬ 
teristic phases occur in different galaxies at different times. 
The gas in the inner region is mostly cold. The cusp forms 
in a wet compaction event that occurs preferentially at in¬ 
termediate redshifts (2 ~ 2 — 4) within the central 1 kpc. 
As mentioned in the Introduction, minor mergers, counter¬ 
rotating streams, recycling and tidal compression drive the 
gas compaction, possibly through VDI, into a compact, star¬ 
forming blue nugget (Dekel & Burkert 2014; Zolotov et al. 
2015). The SFR profile follows the gas profile, in accordance 
with the Kennicutt relation between gas mass and SFR that 
was built-in into the simulations, leading to a central cusp in 
the SFR profile, which characterizes the blue nugget phase. 
In the phase of gas shrinkage and blue nugget formation, 
the extended disc gas mass and SFR is reduced or stays 
constant, while the gas density and SFR in the core goes 
up, indicating an early outside-in quenching phase. The stel¬ 


lar mass density within 1 kpc grows to a maximum in this 
stage, and it remains constant thereafter, i.e., the central 
stellar mass density saturates. As highlighted by Z15 (see 
their Figure 1), based on the same sample of simulations 
used here, the fraction of in-situ star formation in the bulge 
is high, gradually declining from above 60% at 2 ~ 5 to 
~ 50% at 2 ~ 1. This is clear evidence for wet compaction 
at high redshifts. 

Following compaction and the blue nugget phase, the 
centres of the galaxies get depleted from gas, a central hole 
in the gas density develops, and most of the gas is found in 
an outer ring of radius ~ 3 — 10 kpc, and in one case (V07) 
even 15 kpc. Therefore, in this stage, the star-formation also 
takes place in the outskirts, leading to stellar mass growth in 
the outer regions. The galaxies decrease the sSFR within the 
central 1 kpc by 2 — 3 orders of magnitude, while the sSFR 
in the outskirts stays at a similar level as at earlier times. 
We interpret this as inside-out quenching because the star- 
formation activity reduces from inside-out. We will analyse 
the inside-out quenching and its connection to the charac¬ 
teristic evolutionary pattern of our high -2 galaxies in more 
detail in Sections 3.4 and 4, respectively. 

For the low-mass galaxies (Figure 1, left panel), the fea¬ 
tures described before are milder. The average surface den¬ 
sity profiles of gas mass show only little evolution between 
2 = 1 — 4. However, the individual low-mass galaxies show 
substantial evolution of the gas profiles with several phases 
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of high central gas density. The cuspy phases are washed 
out in the stacked profile because these phases happen at 
different epochs for different galaxies. The low-mass galaxies 
tend to grow their stellar mass in a self-similar way, increas¬ 
ing the stellar mass at all radii with the same multiplicative 
pace that keeps the profile shape the same in log-log space. 
This is consistent with the SFR profile roughly following the 
shape of the stellar mass profile, which leads to a rather flat 
sSFR profile out to ~ 7 kpc, especially so at « > 2. Since the 
sSFR is the inverse of the mass e-folding time, the flat sSFR 
profile implies that at all radii, the galaxy roughly doubles 
its mass within the same time interval. Only at late cosmic 
times {z ~ 1) do the low-mass galaxies begin to show a ris¬ 
ing sSFR from the central ~ 0.5 kpc to the outskirts (~ 4 
kpc). 

Summarizing, we find in each galaxy and in the stacked 
profiles at least one phase of wet gas compaction with the 
formation of a compact, gas-rich, star-forming nugget char¬ 
acterized by a cuspy profile of gas and SFR. The blue nugget 
phase is the onset of central gas depletion. Thus, the outside- 
in quenching phase is followed by an inside-out quenching 
phase, associated with stellar mass growth in the outskirts 
at late times. This evolutionary pattern is more noticeable 
in the stacked profiles of high-mass galaxies, since they tend 
to go through just one major compaction phase during the 
same cosmological epoch. Lower-mass galaxies tend to go 
through several compaction events, not all at the same time, 
which washes out the signature. In Section 4, we will stack 
all galaxies according to this evolutionary pattern of com¬ 
paction and quenching, and find that the low-mass galaxies 
also follow a similar evolutionary pattern. 

3.2 Evolution of the Sersic Index n 

We fit the stellar mass surface density profiles with Sersic 
(1968) profiles to characterize their shapes. The evolution 
with cosmic time and stellar mass of the median Sersic index 
n is shown in Figure 2. Overall, galaxies at later cosmic 
time (lower redshifts) and higher masses have a higher Sersic 
index: n increases from 1.5 — 2.0 at M* < 10® Mq and 
2 >4ton = 3 — 4at M* > 10^® Mq and z < 2. This 
shows that the simulated galaxies transform from a disc-like 
n ~ 1 profile to a spheroidal-like n ~ 4 profile. The galaxies 
of the high-mass subsample evolve ahead of the galaxies of 
low-mass subsample and have on average a higher median 
Sersic index. Furthermore, the median Sersic index of the 
high-mass galaxies shows a clear peak of n = 4.61®'g at 
2 ~ 4, indicating that these galaxies are very compact. After 
this compaction, the Sersic index decreases to n = S.dlg'j 
at 2 ~ 1. These classical, de-Vaucouleurs-like profiles are 
consistent with the results from previous simulations with 
weaker feedback models (Ceverino et al. 2015a). We also 
fitted a two-component Sersic model to the mass profile, but 
the two-component fits do not represent the data better, i.e., 
we do not find a stellar surface-density excess towards the 
centres beyond the already rather steep Sersic profile. 

3.3 Inner versus Outer Evolution 

In Figure 3, we investigate the stellar mass growth, either 
within a sphere of 1 kpc or for the whole galaxy within a 



Figure 2. Sersic index n versus stellar mass M*. Squares and 
circles show the median Sersic index for the low- and high-mass 
galaxies, respectively. The colours correspond to redshift accord¬ 
ing to the colour bar. The errorbars indicate the 1 a scatter. The 
median Sersic index increases with stellar mass towards lower red- 
shifts from n ^ 1.5 — 2.0 to n ^ 3 — 4. For the high mass galaxies, 
the median Sersic index peaks (n ^ 4.5) at 2 ~ 4.0. 
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Figure 3. Growth curve of stellar mass for the whole galaxy (blue 
line) and within 1 kpc (red line). The top and bottom panels refer 
to the low-mass (M* < 10^®'® Mq at 2 = 2) and high-mass (M* > 
10^®'^ Mq at 2 = 2) subsamples, respectively. The masses are 
normalized to their value at 2 = 1. The shaded region indicates 
the 1 a scatter. The growth at high redshift is exponential in 2, 
as predicted by Equation 1. For high-mass galaxies, the growth 
of the central stellar mass saturates starting at 2 Ri 3. 
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radius of 0.2x_Rvir, for low-mass and high-mass galaxies. The 
growth curves are normalized by the corresponding stellar 
mass at 2 = 1. We find that low-mass galaxies (top panel) 
increase their mass inside both radii at roughly the same 
pace: the central mass grows concurrently with the total 
stellar mass, though there is a tendency for the mass in the 
central 1 kpc to grow somewhat earlier, e.g., by ~ 0.1 dex 
at 2 = 2 — 3. The mass increases by a factor of ~ 8 from 
2 = 3 to 2 = 1. For the high-mass galaxies (bottom panel), 
the mass in the inner 1 kpc is in place signihcantly earlier 
than the total mass, such that at any redshift in the range 
6 — 2.5, the core mass is a higher fraction of its final mass at 
2 = 1, by a factor of ~ 3. The central and total stellar mass 
grow at a similar pace only at 2 > 3. After 2 ~ 3, the central 
stellar mass saturates and it remains roughly constant after 
2 ~ 2.5, while the total stellar mass continues to grow, by a 
factor of ~ 2 from 2 = 2.5 to 2 = 1. 

We quantify the rate of growth where it resembles a 
straight line in this semi-log plot by htting the curve of 
growth with the theory-motivated functional form 

M( 2 ) = Mo X exp[— q :(2 — 2 o)], (1) 

where the mass at some fiducial redshift 2 o is given to be Mq. 
For a given choice of 2 o, we determine the best-fit value of 
the slope a within a redshift range between 2 o and a certain 
higher 2 . We set 2o = 1. 

The functional form in Equation 1 is motivated by the¬ 
ory, see Dekel et al. (2013) for a simple derivation. Briefly, it 
can be shown that the average specific mass accretion rate 
(sMAR) onto haloes can be approximated in the Einstein- 
deSitter regime (2 > 1) by sMAR = M/M ~ s • (1 -I- 2 )^ 
with n —>■ 5/2. This can be simply integrated, resulting in 
the growth of halo mass as a function of 2 as in Equation 1 
with a = (3/2)sti and ti = H//^ ~ 17.5 Gyr for 

standard ACDM. This functional form was indeed found to 
be a good fit to the growth in earlier cosmological A-body 
simulations of lower resolution (Wechsler et al. 2002; Neis- 
tein & Dekel 2008). Then, one can assume that, even if the 
accretion rate into the central galaxy may be lower than the 
sMAR into the halo, the specific accretion rates are compa¬ 
rable. This has been confirmed in simulations (Dekel et al. 
2013). Finally, conservation of mass through the “bathtub” 
model (Bouche et al. 2010; Lilly et al. 2013; Dekel & Man- 
delker 2014) yields a quasi-stable state with sSFR ~ sMAR, 
so Equation 1 applies also for the growth of stellar mass. 

In T15a we have investigated the evolution of the main 
sequence of SFGs in these simulations and we have found 
that s ~ 0.046 Gyr“^ at M* = 10^° Mq, resulting in a ~ 
1.21. The parameter s has a weak mass dependence oc 
associated with the slope of the AGDM power spectrum on 
the relevant scales. This crudely translates to a similar mass 
dependence in a. Wechsler et al. (2002) found for the halo 
mass growth a ~ 1.2 at Mvir = 10 ^^, and a mass dependence 
of a oc As mentioned in Section 2, our galaxies have 

halo masses of Myir = 2 x 10 ^^ — 2 x 10 ^^, i.e., we expect a 
in the range of 0.97 to 1.31. 

For the low-mass galaxies, we ht a in the redshift range 
2 = 1 — 6 while Hxing 20 = 1.0. For the total stellar mass, 
we find a = 0.99±0.01, i.e., low-mass galaxies increase their 
stellar mass by an e-fold over one unit in redshift. For the 
stellar mass within 1 kpc, we hud a slightly slower growth 
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Figure 4. Evolution of the sSFR for the whole galaxy (blue line) 
and within 1 kpc (red line). The top and bottom panels refer 
to the low-mass (M* < 10^®'^ Mq at 2 = 2) and high-mass 
(M* > Mq at 2 = 2) subsamples, respectively. The shaded 

region indicates the 1 a scatter. The sSFR is rather constant down 
to 2 = 3.5 (4.2) for the low-mass (high-mass) sample, after which 
it declines gradually, rather steeply for the massive galaxies at 
2 s; 2. For low-mass galaxies the sSFR decreases with the same 
pace at all radii, while for high-mass galaxies the sSFR decreases 
faster in the inner regions. 

with a = 0.82 ± 0.03. For the high-mass sample, we can 
already see in Figure 3 that the growth slows down with 
cosmic time, i.e., a decreases. Fitting the data in two redshift 
bins (2 = 1.0 — 3.5 and 2 = 3.5 — 6.0), we find for the total 
mass a = 0.98 ± 0.03 with zq = 3.5 and a = 0.41 ± 0.03 
with 2 o = 1.0 for the high -2 and low -2 bin, respectively. For 
the stellar mass within 1 kpc, we find again only a slightly 
smaller a with a = 0.92 ± 0.04 at high redshifts, whereas 
a = 0.12 ± 0.03 at low redshifts. 

Summarizing, we see that low-mass galaxies throughout 
cosmic time and high-mass galaxies at early times (2 > 3.5) 
grow stellar mass with the same pace in the inner (< 1 
kpc) region and overall (with a semi-log slope a ~ 0.8 — 
1.0). This is in good agreement with the halo growth rate 
of A-body simulations and shows that the sSFR is indeed 
tightly coupled to the specific halo mass increase, as seen 
for an earlier suite of the simulations in Dekel et al. (2013). 
The most interesting result for our purpose here is that, 
for massive galaxies, after 2 ~ 3, the growth rate of stellar 
mass slows down compared to the halo growth rate. This 
is especially pronounced inside the inner 1 kpc, where the 
stellar mass saturates after 2 ~ 3, establishing a compact 
stellar core already then, while the stellar mass continues to 
build up outside this core at later times. 

3.4 Inside-Out Quenching 

Figure 4 shows the sSFR in the central 1-kpc region and 
for the whole galaxy as a function of redshift. Low-mass 
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Figure 5. Half-mass size Re versus stellar mass M*. Squares and 
circles show the median size for the low- and high-mass galax¬ 
ies, respectively. The colours correspond to redshift according to 
the colour bar. The errorbars indicate the 1 a scatter. The gray 
solid and dashed lines indicate the half-light sizes for early type 
(mostly quiescent) galaxies and late-type (mostly star-forming) 
galaxies from observations of Shen et al. (2003) at 2 = 0. The 
corresponding black lines show the half-mass sizes by correcting 
the half-light sizes by 25%. The median size of low-mass galaxies 
stays within 1 — 2 kpc, whereas the one of high-mass galaxies de¬ 
creases from 2 = 4.9 to 4.0 from 1.5 kpc to 0.6 kpc and increases 
from 2 = 4.0 to 1.0 from 0.6 kpc to 2.2 kpc. 

galaxies have a constant sSFR down to 2 = 3.4, after which 
it declines gradually. The sSFR is the same within the two 
apertures at 2 > 2.5. At lower redshifts we see a deviation: 
the sSFR is declining, slightly faster in the central 1 kpc than 
over the whole galaxy. The average total sSFR is decreasing 
from 3 Gyr“^ at 2 ~ 6 to 0.3 Gyr“^ at 2 ~ 1, while the 
sSFR within 1 kpc drops from 3 Gyr“^ to below 0.2 Gyr“^. 

High-mass galaxies have a constant sSFR till 2 = 4.2, 
before the sSFR declines. Below 2 = 3.5, the sSFR declines 
faster in the inner 1 kpc region. The sSFR in the central 
region is 3.2 Gyr“^ at 2 ~ 6 and it decreases to 0.01 Gyr“^ 
at 2 ~ 1, while the total sSFR of the whole galaxy only 
decreases to 0.07 Gyr”'^ at 2 ~ 1. The sSFR reduces in the 
centre by more than two orders of magnitudes, while the 
sSFR in the outskirts stays high and declines later. This 
is consistent with a picture of inside-out quenching in the 
post-compaction phase and it is in line with what we have 
seen in Figure 3, where the stellar mass within the central 
1 kpc saturates after 2 ~ 3, while the mass in the outskirts 
keeps growing. 

3.5 Size-Mass Evolution 

Z15 studied the size evolution of individual galaxies in the 
simulations presented here (see their Figure 9). They found 
that the galaxies indeed go through a compact, nugget 
phase, with a compactness that is comparable to observa¬ 
tions. We take this discussion one step further and look at 


the median effective size of all galaxies in our evolving sam¬ 
ple, as a function of stellar mass and redshift, and compare 
it with observations. An important point is that most sizes 
quoted in observations refer to the half-light radius, whereas 
we determine the stellar half-mass radius in the simulations 
(and do it in three dimensions). Szomoru et al. (2013) (see 
also Hernquist 1990) found that the average (rest-frame g- 
band) light-mass size difference of ~ 25% is roughly the same 
at all redshifts (0.5 < 2 < 2.5), and does not strongly corre¬ 
late with stellar mass, specific star formation rate, effective 
surface density, Sersic index, or galaxy size. 

In Figure 5, we plot the median half-mass size Re versus 
stellar mass M* for our low-mass (squares) and high-mass 
(circles) samples in the redshift range 2 = 1 — 6. At high 
redshifts (2 > 4) and low stellar masses (M* < 10® Mq), 
both low-mass and high-mass galaxies increase their sizes 
according to 

AlogRe = 0.3AlogM*, (2) 

i.e., galaxies increase their size by a factor of 2 for a factor 
of 10 increase in their stellar mass. At 2 ~ 3 — 4, high-mass 
galaxies show a significant contraction in a relatively nar¬ 
row mass range around M* ~ 10®'® Mq. The median size 
decrease from 2 = 4.9 to 4.0 from 1.2'Aq \ kpc to O.SGIq q® 
kpc, and increase from 2 = 4.0 to 1.0 from 0.561 q'q 2 kpc to 
3.0 ('Iq' 3 kpc. The lower mass galaxies also show the signature 
of compaction, but to a lesser extent: they compact by a fac¬ 
tor of ~ 20% near M* ~ 10®'® Mq and 2 ~ 2 — 3. The more 
mild compaction seen on average for the low-mass subsam¬ 
ple indicates that each compaction event is less pronounced 
(e.g., because it happens later when the density is lower and 
the gas fraction is lower), and that the compaction events 
are spread over a larger range of redshifts and masses, so 
their signal is smeared in the stacked evolution track. 

In the post compaction phase, the galaxies increase 
their size more significantly. Their growth in the mass- 
size plane is well approximated by A log Re = 0.5AlogM* 
and A log Re = O.SAlogM* for the low-mass and high- 
mass galaxies, respectively. Both slopes are clearly steeper 
than before compaction. Due to the pronounced compaction 
phase at intermediate redshifts and masses, there is only a 
mild overall growth of the median size from 2 = 6 to 1, 
where Re grows by only a factor ~ 4 as the mass grows by a 
factor ~ 100. The overall growth of size is minor compared 
to the mass growth for all galaxies over the same period. 
This is consistent with the picture discussed before in which 
the stellar mass increases at all radii with a similar pace. 
Furthermore, Re for the more massive galaxies tends to be 
smaller at every given mass, indicating that galaxies that 
formed earlier in a denser universe have also on average a 
higher stellar (effective) density (e.g., Carollo et al. 2013a). 

Interestingly, the maximum compaction tends to oc¬ 
cur in a relatively small range of stellar masses around 
M* ~ 10®'® Mq. This is apparently the stellar mass range 
where the core switches from being dark-matter-dominated 
to baryon-dominated, which largely occurs near the point of 
maximum compaction (Zolotov et al. 2015; Geverino et al. 
2015b; Tomassetti et al., in prep.). 

In Figure 5, the solid and dashed lines indicate the crude 
observed size evolution of early-type and late-type galaxies 
at 2 = 0 of Shen et al. (2003). The gray lines show the 
measured half-light sizes, whereas the black lines show the 
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Figure 6. Average fraction of the stellar mass in the outskirts 
(1 — 5 Re) formed in-situ and ex-situ. A large fraction of ~ 80% of 
the stellar mass in the disc has formed in-situ, i.e., in the galaxy, 
and only about 20% has been accreted. This fraction stays roughly 
constant with redshift. 


half-mass sizes (assuming a constant correction factor of 25 
%). We see that the sizes of our low and high-mass galaxies 
at z = 1 are consistent with the sizes of early-type galaxies 
at z = 0, which are mostly quiescent galaxies. As shown 
before, most of our simulated galaxies are in process of being 
quenched at z = 1. Therefore, the simple shut down of star- 
formation in our galaxies with no additional size (and mass) 
growth would roughly produce the observed galaxies at z = 
0 . 

The late {z < 3) size growth of our simulated galaxies 
can be understood by the addition of stellar mass in the out¬ 
skirts, as shown in Section 3 by analysing the stellar density 
profiles. To understand if this outer stellar mass growth hap¬ 
pens due to accretion of stellar mass or in-situ star formation 
in an extended ring, we plot in Figure 6 the fraction of stel¬ 
lar mass of the stars in the outskirts that formed in-situ 
versus ex-situ. At a certain redshift, stars in the outskirts 
were selected by choosing all stars within 1 — 5 7?e. If a star 
was formed outside of 10% of the virial radius (> 0.1 i?vir), 
then the star was classified as ex-situ. On the other hand, all 
stars that formed within 0.1 i?vir were classified as in-situ. 
We find that most stellar mass in the outskirts formed in- 
situ: ~ 80% of the mass was formed in-situ at all redshifts 
down to 2 « 1. 


4 EVOLUTIONARY PHASES 

We are now going to interpret the evolution of the sur¬ 
face density profiles in the context of the characteristic evo¬ 
lutionary phases of high-z galaxies. As mentioned in the 
Introduction, Z15 put forward a picture in which galax¬ 
ies evolve through gas compaction to compact, star-forming 
blue nuggets and then quench inside-out. Instead of stacking 
the profiles at a certain redshift as carried out in Section 3, 
we now stack the profiles of the galaxies at four evolution¬ 
ary stages: the pre-blue nugget phase, i.e., when the galaxy 
was an SFG close to the MS ridge, the blue nugget phase 
of peak central gas density and SFR, the post-blue nugget 
phase when central quenching has started, and the quench¬ 


ing phase, where the core has quenched to a red nugget. The 
individual profiles of the simulations are stacked to have the 
same effective density, following the procedure outlined in 
Section 3.1. Figure 7 shows from top to bottom the surface 
density profiles of stellar mass, gas mass, SFR and sSFR. In 
the left panels, we plot the profiles as a function of the phys¬ 
ical radius in kpc, while in the right, the profiles are shown 
as a function of the radius with respect to Re- The differ¬ 
ent lines indicate the mean profile at the pre-blue nugget 
phase (“pre-BN”), the blue nugget phase (“BN”), the post¬ 
blue nugget phase (“post-BN”), and the quenching phase 
(“quench”). 

These four phases have been identified for each galaxy 
individually based on Z15 and T15a. Briefly, in T15a we 
studied in detail the galaxy properties with respect to the 
star-forming main sequence ridge. The blue nugget phase is 
the snapshot in which the galaxy reaches its global maxi¬ 
mum in the central 1 kpc gas mass. This happens for all 
galaxies at 2 : ~ 1 — 4. The pre-blue nugget phase is iden¬ 
tified to be the snapshot where the galaxy is still on the 
main sequence ridge (distance from the main sequence <0.1 
dex), before leaving it towards higher sSFR and entering the 
blue nugget phase. This is typically prior or during the com¬ 
paction process. Usually, this is 3 — 5 snapshots (300 — 500 
Myr) before the blue nugget phase. The post-blue nugget 
phase is defined to be the snapshot after the blue nugget 
phase, where the galaxy is arriving at the lower edge of the 
main sequence, which usually is 2 — 4 snapshots (200 — 400 
Myr) after the blue nugget phase. Finally, the quenching 
phase is the point in time where the galaxy has the lowest 
sSFR with respect to the main-sequence ridge, which usually 
is one of the final snapshots (2 ~ 1 — 2). 

In Figure 7, there is a clear gas cusp in the inner 1 
kpc in the blue nugget phase, where the gas surface density 
reaches > 10® Mq kpc“^ in the centre. Before and after the 
blue nugget phase (pre-BN and post-BN), the gas surface 
density in the inner region amounts to ~ 10^ Mq kpc“^. 
Therefore, a gas compaction event takes place where the 
inflow rate > SFR -f outflow rate, i.e., the contraction is 
wet. The SFR profile follows closely the gas mass profile, 
following the SFR recipe incorporated in the simulations. 
During the blue nugget phase, the SFR density in the centre 
increased by an order of magnitude and the SFR profile is 
cuspy. A weak, but interesting feature is that the sSFR in the 
disc before compaction is high (> 0.01 Gyr“^) out to 20 kpc, 
whereas after the blue nugget phase, this high sSFR regime 
in the disc has shrunk to about 10 kpc. This is explained 
by the fact that in the compaction phase, the SFR slows 
down in the extended disc, because its gas has contracted 
and not replenished yet, while in the centre the SFR goes 
up to a peak at the blue nugget phase. We interpret this 
as outside-in quenching in the early, pre-blue nugget phase. 
After the blue nugget phase, the flow rate in the central 1 
kpc tips to inflow rate < SFR -|- outflow rate, reducing the 
central gas mass and SFR density. In the pre-BN, BN, and 
post-BN phases, the sSFR profile is flat out to ~ 7 kpc, and 
therefore, the surface stellar mass density profiles are nearly 
identical in log-log shape, showing again self-similar growth. 

During the quenching phase, the gas in the central re¬ 
gion gets depleted and a gas hole, usually associated with an 
outer ring, develops. This is consistent with the picture of 
Z15 (see their Figure 6), where a post-compaction extended 
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Figure 7. Stacked profiles during the four galaxy evolution phases. Shown from top to bottom are the stellar mass density, gas mass 
density, SFR density, and sSFR. On the left, the radius is in kpc, and on the right it is scaled by the effective radius Re. The individual 
profiles of the simulations are stacked to have the same median effective density. The four phases are pre-blue nugget (pre-BN; dashed 
purple), blue nugget (BN; solid blue), post-blue nugget (post-BN; dashed green), and quenching phase (quench; solid red), where the red 
nugget is surrounded by a star-forming ring. The shaded areas show Icr scatter. 
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gas ring forms, in which the SFR resumes. This is proba¬ 
bly from fresh gas that comes with high angular momentum 
and cannot get rid of it. At later times, the SFR density 
varies only a little in the outskirts, whereas in the centre, 
the SFR density declines by about two orders of magnitude. 
The sSFR shows the typical signature of inside-out quench¬ 
ing, where the sSFR in the central part is substantially lower 
(by about 2 orders of magnitude) than in the outskirts. The 
absolute difference in the stellar mass surface density profile 
is small, highlighting the convergence of the mass profile at 
later evolutionary time, after the compaction and through¬ 
out the quenching process. 

Comparing the profiles as a function of radius in kpc 
with the radius scaled by Re, there are no qualitative differ¬ 
ences. We see that the cusp of gas density and star-formation 
rate develops inside a radius of r < 2 kpc, which corresponds 
to < 2 Re. The sSFR profiles rises in the red-nugget phase 
out to ~ 5 kpc, which corresponds to about 2 Re. 

Summarizing, a wet compaction causes a development 
of a gas central cusp in the blue nugget phase. There is an 
associated peak of SFR in the centre, while the SFR at larger 
radii decreases (outside-in quenching). The stellar mass pro¬ 
file grows self-similarly, consistent with the flat sSFR profile. 
After the blue nugget phase, central gas depletion leads to 
central quenching of the SFR, leading to the saturation of 
the stellar core. The sSFR profile is rising in the post-blue 
nugget phases, which can be termed inside-out quenching. 


5 COMPARISONS TO OBSERVATIONS 

Thanks to the Hubble Space Telescope {HST), accurate 
measurements of the stellar mass density on scales of ~ 1 
kpc have become possible. Several studies have investigated 
the stellar mass density distribution within galaxies at high 
redshifts (e.g., Elmegreen et al. 2009; Forster Schreiber et al. 
2011; Wuyts et al. 2012; van Dokkum et al. 2013; Patel et al. 
2013; Lang et al. 2014; Morishita et al. 2015). In addition to 
stellar mass density profiles, T15b presented resolved star- 
formation rate density distributions in 2 ~ 2.2 galaxies. In 
this section, we compare our simulated stellar mass and SFR 
density profiles to these observations. 

5.1 Central Stellar Mass Density 

In Figure 8, we plot the central stellar mass density within 1 
kpc, Em*,i, as a function of stellar mass, M*, in our simula¬ 
tions compared with observations. The colour-coding of the 
symbols in the left panel corresponds to redshift. We find a 
rather tight, nearly linear relation between Em*,i and M*. 
All galaxies grow along a well-defined sequence during their 
SFG phase, i.e., they are increasing the central stellar mass 
density concurrently with the total stellar mass, following 
the relation Ea-/*,! oc jyf i o 2 ±o. 05^ Ceverino et al. (2015a) 
used different simulations with much weaker feedback, but 
finding the same linear relation between Em*,i and M*, i.e., 
this suggests that the overall trend is independent of feed¬ 
back. 

At the high-mass end, around ~ Mq, the evolu¬ 

tion tracks tend to flatten: the central stellar mass density 
saturates, while the total mass continues to increase (see also 
Fig. 16 of Z15). We find the relation Eai*,! oc for 


galaxies with M* > Mq. This flattening starts at 

the end of the compaction phase and the tracks remain fiat 
during the quenching process. This implies that any further 
stellar mass growth happens in the outskirts, avoiding the 
cores, consistent with the aforementioned picture. 

The colour-coding of the symbols in the right panel of 
Figure 8 corresponds to the total sSFR. We compare our 
simulations at 2 = 1 — 4 with observations at 2 ~ 2 and 
at 2 ~ 0. As the 2 ~ 2 comparison sample, we use the 
zC-SINF sample of Tacchella et al. (2015c), indicated by 
the star symbols. At 2 ~ 0, we show the best-fit line to the 
quiescent sample of SDSS from Fang et al. (2013) as red solid 
line (Eat,,! oc The contours show the 2 ~ 0 ZENS 

sample (Carollo et al. 2013b; Cibinel et al. 2013a,b) of star¬ 
forming (blue) and quiescent (red) galaxies. The zC-SINF 
sample as well as our simulations are not mass complete, 
and therefore, the comparison has to be take with care. 

We find excellent agreement between the overall shape 
of the Ea/*,! — M* sequence in the simulations and in the ob¬ 
servations, both for SFGs and for quenched galaxies. T15b 
already pointed out that galaxies at high and low redshifts 
practically lie on top of each other in the Ea^*,! — M* di¬ 
agram. This may indicate that galaxies grow along this se¬ 
quence, which is indeed confirmed in the simulations (Fig¬ 
ure 8 and also Z15). In both the simulations and observa¬ 
tions, there is a flattening of the Evr, ,1 — M* relation in 
the high-M*/high-E at*,! end, associated with the central 
quenching process. 


5.2 Density Profiles 

The total stellar masses and SFR are two fundamental global 
parameters of each galaxy and therefore have been measured 
in many observational works. However, only very recently, 
it became possible to resolve galaxies with broad-band pho¬ 
tometry from the HST at high redshifts (2 ~ 2 — 3) and 
determine stellar masses on spacial scales of 1 — 2 kpc (e.g., 
Wuyts et al. 2012; van Dokkum et al. 2013; Patel et al. 2013; 
Lang et al. 2014; Morishita et al. 2015). Morishita et al. 
(2015) used HST /WFG3 and ACS multi-band imaging data 
taken in CANDELS and 3D-HST to study the stellar mass 
growth of the 2 = 0.5 — 3 progenitors of local galaxies with 
mass comparable to the Milky-Way and more massive galax¬ 
ies, based on a constant cumulative number density selection 
of galaxies ^. They conducted a radially resolved pixel SED 
fitting to obtain the radial distributions of the stellar mass 
and rest-frame colours by stacking several tens of galaxies 
at a given redshift and mass bin. They found that the stel¬ 
lar mass profiles of Milky-Way-like progenitors increase in 
a self-similar way, irrespective of radial distance. However, 
the profiles of massive galaxies grow inside-out; about 75% 
of their current total mass has grown since 2 ~ 2.5 outside 
the inner 2.5 kpc. This is consistent with previous studies 


^ Connecting progenitors and descendants using a fixed comov¬ 
ing number density arguments is problematic, especially over long 
time-scales such as form 2 = 2 and 2 = 0. More careful analyses 
using an evolving number density approach or numerical simu¬ 
lations yield 2 — 3 times lower masses for the progenitors (e.g., 
Behroozi et al. 2013b; Lu et al. 2015; Torrey et al. 2015). 
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Figure 8. Central stellar mass density within 1 kpc (Sm*,i) a function of stellar mass (M*). The grey shaded area above the dashed 
line marks the forbidden parameter space where ^ M*/ 7 r. The colour coding in the left panel corresponds to redshift, while the 

colour coding in the right panel indicates the total sSFR, highlighting the difference between SFGs and quiescent galaxies. All galaxies 
grow their central stellar mass density hand-in-hand with their total stellar mass. At M* > Af 0 , there is a flattening in the 

— Af* relation, where the total stellar mas still increases while the central stellar mass density has saturated. We compare our 
simulated galaxies with observations at z ~ 0 and 2 ~ 2 in the right panel. In particular, the blue and red contours show the distribution 
of 2 ~ 0 star-forming and quiescent galaxies of ZENS (Carollo et al. 2013b). The red solid line marks the best fit to the 2 ~ 0 red 
sequence galaxies of SDSS (Fang et al. 2013). As the 2 ~ 2 comparison sample, we use the zC-SINF sample of T15b, plotted with star 
symbols. There is a good agreement between observations and simulations, and between 2 ~ 2 and 2 ~ 0. 


conducted by van Dokkum et al. (2013) and Patel et al. 
(2013). 

As shown in Section 3.1, we find exactly the same qual¬ 
itative behaviour in our simulations: the low-mass galaxies 
preserve the log-log shape of the stellar mass density profile 
at all redshifts, whereas the more massive galaxies grow self¬ 
similar ly only at 2 > 3. Toward lower redshifts, the stellar 
mass growth saturates first in the centre. The stellar mass 
in the outskirts, and thus the total stellar mass, still in¬ 
creases at lower redshifts; ~ 50% of the total stellar mass 
is assembled between 2 ~ 2 and 2 ~ 1 in the outskirts 
(r > 2.5 kpc). More quantitatively, we compare the expo¬ 
nential growth rate a from Equation 1 as derived from the 
observations and the simulations. For the lower mass galax¬ 
ies, we fit the data at 2 = 1.0 — 2.5 with 20 = 1.0. We find 
a = 1.16 and a = 1.17 for the data of van Dokkum et al. 
(2013) and Morishita et al. (2015), respectively. In the sim¬ 
ulations (Section 3.3), we obtained a = 0.99, which is in the 
same ball park, perhaps slightly smaller, i.e., the lower mass 
galaxies in the simulations tend to form stars a little ear¬ 
lier then the ones in the observations. For the more massive 
galaxies, we fit the data of Patel et al. (2013) at 2 = 1.0 —2.5 
with 2o = 1.0, finding a = 0.52. The data of Morishita et al. 
(2015) shows a break at 2 = 2.0, so we fit two redshift bins 
independently from each other. At 2 = 2 — 3 with 20 = 2 


we find a = 1.89, and for 2 = 1 — 2 with 20 = 1 we find 
a = 0.46. In the simulations, recall that we find a = 0.41 
with 2 o = 1.0 at 2 = 1.0 — 3.5, which is also in reasonable 
agreement with the observations. 

Morishita et al. (2015) used the rest-frame U — V colour 
to infer that the quenching of the central dense region is 
ahead of the outer disc, which is also qualitatively consis¬ 
tent with the inside-out quenching found in the simulations. 
Unfortunately, the rest-frame U — V colour is an uncer¬ 
tain tracer of SFR, since it depends on the assumed star- 
formation history as well as dust properties. Instead of us¬ 
ing broad-band photometry as a tracer for star-formation, 
T15b used the Ha emission line from VXT/SINFONI inte¬ 
gral field unit spectroscopy as part of the SINS/zC-SINF 
effort (Newman et al. 2013; Forster Schreiber et al. 2014; 
Genzel et al. 2014; Tacchella et al. 2015c, b). The Ha emis¬ 
sion line traces the recent SFR (on timescales of 10^ yr), 
independent of the previous star-formation history. Thanks 
to adaptive optics, the same HST-like, 1 kpc resolution is 
obtained for the Ha emission line maps. T15b presented a 
sample of 22 2 ~ 2.2 galaxies with both HST stellar masses 
and SINFONI SFR maps, which spans a wide range in stellar 
mass (4 x 10® — 5 x 10^^ Al©) and SFR (20 — 300 M© yr“^), 
probing the star-forming main sequence at 2 ~ 2. 

In Figure 9 we compare the stacked stellar mass, gas 
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Figure 9. Comparison of the shapes of the stacked profiles from 
the simulations at z = 1.5 — 3.0 with the estimates from ob¬ 
servations of star-forming galaxies (Tacchella et al. 2015b). The 
stacked surface density profiles of stellar mass, gas mass, SFR, 
and sSFR are shown from top to bottom. The color coding corre¬ 
sponds to three different mass bins as indicated in the legend. The 
circles, connected by dotted lines, show the observational data. 
The thick solid lines show the renormalized simulated galaxies, 
whereas the dashed lines show the measured profiles of the sim¬ 
ulated galaxies. Overall, we find good agreement in the profile’s 
shapes between simulations and observations. However, the nor¬ 
malization, especially of the SFR profiles, differ up to ~ 1 dex 
between observations and simulations. 


mass, SFR and sSFR profiles of the simulations with the 
observations at 2 1.5 — 3.0. We stack the simulated galax¬ 
ies in the same mass bins as the observations of T15b 
(with a median redshift of 2 ~ 2 . 2 ), namely, in a low- 
(9.3 < logj^Q M*/Mq < 10.1), an intermediate- (10.1 < 
logj^Q M*/M 0 < 10 . 6 ), and a high-mass bin ( 10.6 < 
logj^Q Mi,/MQ < 11.5). The individual profiles of the simu¬ 
lations and observations are each stacked to have the same 
effective stellar mass (and SFR) density, following the de¬ 
scription outlined in Section 3.1. The dashed lines show the 
stacked profiles of the simulations. In order to compare the 
profiles’ shapes, the main focus of this comparison presented 
here, we may need to re-normalize the profiles of the simu¬ 
lations by requiring them to have the same Eefi as the one 
from the observations (shown as solid lines). 

Important to mention are the following two caveats in 
the comparison between simulations and observations in Fig¬ 
ure 9. First, the SFRs are estimated in different ways, both 
associated with significant uncertainties. While the observed 
SFR is estimated from the dust-corrected Hq emission-line 
flux, in the simulations it is estimated by counting the newly 
born stars over a given period. Second, both the observed 
and simulated samples are not complete samples in terms of 
mass or SFR, incomplete in different ways, which could lead 
to a bias. In particular, the observed sample consists mainly 
of SFGs on the main sequence, while no explicit selection 
on SFR is applied in the simulations. This can explain at 
least part of the difference between the average SFR at a 
given mass, which is lower in the simulations compared to 
the observations. The motivation for the comparison of the 
simulations to observed SFGs is our finding in T15a (where 
we modeled the star-forming main sequence self-consistently 
within the simulations) that most of our simulated galaxies 
are on the main sequence at 2 > 2. However, they do tend 
to fall bellow the Main Sequence toward lower redshifts and 
when they have higher masses, thus introducing a bias in 
the comparison. What matters to us here is that the shapes 
of the profiles in the different phases of evolution, including 
the SFR profiles, are consistent between the simulations and 
observations. This indicates that the simulations capture the 
important physical processes, at least qualitatively, and al¬ 
low us to understand the evolution through compaction and 
quenching phases. 

The top panel of Figure 9 shows the simulated and ob¬ 
served stacked stellar mass density profiles. We obtain as a 
re-normalisation factor for the stellar mass profiles of the 
low-, intermediate-, and high-mass sample small offsets of 
-1-0.13, —0.05 and -fO.lO dex, respectively, namely, the re¬ 
quired re-normalization is negligible. In the intermediate- 
and high-mass bin, we find good agreement in the log-log 
profile shape. In the lower mass bin, there is a certain dif¬ 
ference in the profile shapes: the profile from observations 
shows a more disc-like profile (with n ~ 1 ), whereas the one 
from simulations is more spheroidal-like (with n ~ 3, see 
Fig. 2). 

The low-mass profile of the simulations does resemble 
the high-mass profile of the observations. The second panel 
from the top shows the gas mass profile of the simulations 
only, since there is no available data for the resolved gas 
mass profile. We see a rather flat profile between 1 and 5 
kpc, and a steeper decline at larger radii. 

In the bottom two panels of Figure 9, we show the SFR 
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Figure 10. Progression of the “quenching wave”. The red 
line shows progression of quenching within the massive galaxies 
(M* > 10iO'2 Mq ) in our simulation (shaded regions indicates 
the 1 a scatter). We find that the galaxies quench their outskirts 
(> 5 kpc) 1 — 2 Gyrs after they started quenching in their centre. 
The blue line represents the empirical model based on observa¬ 
tions of Tacchella et al. (2015b), showing excellent agreement. 

and sSFR profiles. We find a substantial re-normalisation 
factor for the SFR profiles: we have to re-normalize the sim¬ 
ulated SFR profiles of the low-, intermediate-, and high-mass 
sample by a factor of -1-1.39 dex, -1-0.91 dex, and +1.33 dex, 
respectively. This offset between observations and simula¬ 
tions could partly be due to the caveats mentioned above, 
and possibly also due to too low gas fractions and total 
SFRs in the simulations at 2 < 2 as a result of the SFR be¬ 
ing overpredicted at much higher redshifts (see also T15a), 
but possibly also due to overestimation in the observations. 

Focusing on the shape of the profiles, we find overall 
good agreement between the profiles’ shapes in the obser¬ 
vations and in the simulations. For the low-mass galaxies, 
the simulated SFR profile is slightly flatter in the inner 3 
kpc than the observed one. In the intermediate-mass bin it 
is the opposite. However, the shapes of the SFR prohles of 
the simulated galaxies agree within the errors with the ob¬ 
served ones. The high-mass prohle of the simulations and 
observations are in very good agreement. As highlighted in 
Section 3, the sSFR is increasing towards the outskirts. This 
effect of reduced star-formation activity toward the centre 
is stronger for intermediate-mass galaxies than for low-mass 
galaxies. As discussed in Sections 3 and 4, the rise in the 
sSFR prohle can be interpreted as a central gas depletion 
after the blue nugget phase, which is followed by an ex¬ 
tended star-forming ring. This is consistent with the inside- 
out quenching scenario proposed by T15b and Z15. Overall, 
the agreement between the shapes of the observed and sim¬ 
ulated prohles at 2 ~ 2 is an encouraging sanity check. It 
indicates that the observed prohle shapes can be interpreted 
as rehecting the build-up of a compact bulge by wet com¬ 
paction and the subsequent inside-out quenching, as seen in 
the evolution of the simulated prohles in Figure 7. 

5.3 Inside-Out Quenching 

To quantify the inside-out quenching phase more quantita¬ 
tively in the simulations and allow a more direct comparison 
to observations, we measure the time when the sSFR falls 
below a critical sSFR threshold at every given radius. We 


set this threshold to be sSFRthresh = 0.1 Gyr“^. The moti¬ 
vation for this number is that a sSFR of 0.1 Gyr“^ roughly 
corresponds to a mass doubling time of 10 Gyr, i.e., several 
times the Hubble time at z >2. This means that the mass 
increase from star-formation below sSFRthresh is negligible. 
Changing sSFRthresh in the range 0.05 — 0.2 Gyr“^ does not 
change our results significantly. 

Figure 10 shows the progression of the inside-out 
“quenching wave”. For all massive galaxies (M* > 

10^°'^ Mq), we measure at which cosmic epoch and ra¬ 
dius the sSFR drops below sSFRthresh- The “quenching” 
time tquench of a galaxy starts ticking when the first radial 
bin quenches. The massive galaxies in the simulations start 
quenching (tquench = 0) at 2 = 2.4to'l. By 2 = 1, seven 
out of 12 massive galaxies are fully quiescent, i.e., they are 
> 0.3 dex below the main sequence (T15a), while the other 
5 galaxies lie on the lower edge of the main sequence. We 
see that in the simulations, the galaxies quench their inner 
regions first, while the outskirts quench 1 — 2 Gyrs later (red 
line in Figure 10). 

In T15b we have used an empirical model to constrain 
the progression of inside-out quenching motivated by obser¬ 
vations. The main idea was to use the measured stellar mass 
and SFR profiles for the most massive galaxies in that sam¬ 
ple (M* ~ 10^^ Mq) and compare them with local 2 ~ 0 
galaxies. By assuming that the galaxies’ star formation rates 
follow the ones of the star-forming main sequence at any 
given time, we estimated when the galaxy has to stop form¬ 
ing stars at the certain radius so that they do not overshoot 
the 2 = 0 stellar mass profile. We found in T15b that 2 ~ 2 
galaxies have to cease the star-formation in their centres 
very soon (<< 1 Gyr), while they keep forming stars for a 
few Gyr in their outskirts. Since this empirical model does 
not take into account mergers or stellar migration, the ob¬ 
tained timescale for quenching is an upper bound, since dry 
merging would contribute to the mass profile as well, with¬ 
out increasing the SFR there. 

In Figure 10 we show the result from this empirical 
model in blue. As one can see, the model agrees with the 
simulations out to a few (~ 5 kpc). At r ~ 10 kpc, the 
simulations predict a shallower slope in tquench than what we 
have predicated with the T15b empirical model. As men¬ 
tioned above, this difference can possibly be explained by 
merging. 

5.4 Physics of Quenching 

What is the physical origin of the inside-out quenching? This 
is discussed in detail in our companion papers (T15a and 
Z15). Specifically, in T15a, we study in detail the flow rates 
in the central 5 kpc (see Figure 9 and D1 of that study). We 
find that the blue nugget phase comes rather quickly to an 
end when the inflow rate into the central regions becomes 
smaller than the SFR + outflow rate inside that region. The 
gas depletes in the centre quickly due to the formation of 
new stars and the associated stellar feedback that derives 
winds, while there is no gas supply from the disc that has 
shrunk and at least temporarily disappeared. During this 
phase, the galaxies reduce its star-formation substantially 
(especially in the inner regions) and the sSFR profile is rising 
toward the outskirts as a new, extended gas ring is forming. 
Newly accreted gas through the halo and into the galaxy 
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can initiate new star-formation, such that the galaxy does 
not necessarily quench for good. 

For full and permanent quenching, we argue in T15a 
that the disc replenishment time (i.e., timescale for accre¬ 
tion of fresh gas into the galaxy) is longer than the depletion 
time (sometimes also called the gas consumption time). This 
is the case when the circumgalactic medium in the host halo 
becomes hot, typically when the halo mass is above the crit¬ 
ical mass for virial shock heating, and at late times when 
cold streams do not penetrate efficiently (Birnboim & Dekel 
2003; Keres et al. 2005; Dekel & Birnboim 2006; Zolotov 
et al. 2015). As mentioned in Section 2.2, full quenching to 
very low sSFR is not always reached in the current simula¬ 
tions because the feedback may still be underestimated (e.g., 
the adiabatic phase of supernova feedback is not resolved, 
and AGN feedback is not yet incorporated). However, 12 
out of our 26 galaxies drop to more than Icr below the main 
sequence ridge hy z = 1, and all are continuously reducing 
their sSFR for the last several Gyrs, indicating that they are 
in a long-term quenching process. 


6 CONCLUSION 

We have used a suite of 26 zoom-in cosmological simula¬ 
tions to investigate the development of the density profiles 
in galaxies at 2 = 1 — 7 through their characteristic phases 
of evolution. The galaxies compactify to blue nuggets and 
then quench inside-out, developing a red-nugget core that 
is usually surrounded by a star-forming ring. Overall, we 
hnd good agreement between observations and simulations 
for the main features of wet compaction, inside-out quench¬ 
ing, and saturation of dense stellar cores. We summarize our 
main findings below. 

To shed light on the physical process in these high -2 
galaxies, we have stacked the proHles of the galaxies based 
on their four evolutionary states with respect to compaction 
and quenching: the pre-blue nugget, blue nugget, post-blue 
nugget, and quenching phases, following Z15 and T15a. The 
emerging picture is that a wet compaction (Dekel & Burkert 
2014), caused by gas-rich mergers or smoother gas streams 
(such as counter-rotating streams and low-angular momen¬ 
tum recycled gas), leads to a high central gas density above 
10® M© kpc“^ and an episode of high central SFR, i.e. a blue 
nugget. This short phase can be interpreted as an outside-in 
quenching episode, where the SFR slows down in the shrink¬ 
ing disc, while in the centre the SFR peaks in the blue nugget 
phase. In this phase, the gas and SFR profiles are cuspy. 

The blue nugget phase comes rather quickly to an end 
when the balance of the flow rates in the central 1 kpc tips to 
inflow rate < SFR -|- outflow rate. The immediate reason for 
this is that while the SFR -f outflow rate are at their max¬ 
imum, the inflow is suppressed because the disc has largely 
disappeared (shrunk). In parallel, the bulge has grown, so 
the remaining disc is stable gravitationally (morphological 
quenching, see Martig et al. 2009; Cacciato et al. 2012; Gen- 
zel et al. 2014; Tacchella et al. 2015b) and does not drive 
more gas inward. After the gas has depleted in the centre, 
a star-forming, gas-rich ring forms in some of the galaxies. 
This is consistent with the ID axisymmetric thin disc model 
of Forbes et al. (2014) based on continuous disc instability 
with Toomre (1964) Q-parameter of the order of unity. They 


showed that at every radius galaxies tend to be in a slowly 
evolving equilibrium state (as predicted in a cosmological 
setting by Dekel et al. 2009b and in Gacciato et al. 2012) 
and they explicitly predicted the formation of a central gas 
“hole” with an extended gas ring. 

In the post-blue nugget phase, the sSFR is rising to¬ 
ward the outskirts, which we interpret as inside-out quench¬ 
ing phase. We argue that a galaxy departs from the main 
sequence and quenches once the disc replenishment time is 
longer than the core depletion time. This naturally hap¬ 
pens at late redshifts, and earlier for more massive galaxies 
(“downsizing”, T15a). Quenching can be maintained for long 
term when the replenishment time becomes even longer, e.g., 
once the circum galactic medium in the host halo becomes 
hot, typically when the halo mass is above the critical mass 
for virial shock heating, and at late times when cold streams 
do not penetrate efficiently (Birnboim & Dekel 2003; Keres 
et al. 2005; Dekel & Birnboim 2006; Zolotov et al. 2015). 

The surface stellar mass density profiles of low-mass 
galaxies grow self-similarly at 2 = 1 — 7; on average, the 
stellar mass density increases at all radii with a similar mul¬ 
tiplicative pace. This can also be characterised by a flat aver¬ 
age sSFR profile. The rate of (total) stellar mass growth for 
the low-mass galaxies is well described by M{z) = M{z = 
1) X exp[— 1 . 0(2 — 1)], which is in excellent agreement with 
the halo growth rate, as measured in simulations and pre¬ 
dicted by simple theory, and shows that the sSFR is tightly 
coupled to the specific halo mass increase (Wechsler et al. 
2002; Neistein & Dekel 2008; Dekel et al. 2013). 

The stellar mass profiles of high-mass galaxies grow self- 
similarly only at high redshifts (2 > 3.5). At lower redshifts, 
the high-mass galaxies decouple from the halo mass accre¬ 
tion and grow signihcantly slower. Comparing the growth 
of the central region with the outskirts, we find that most 
of mass increase after 2 < 3 happens in the outskirts: the 
stellar mass within 1 kpc saturates starting at 2 « 2.5. This 
is consistent with the strong decline with time of the sSFR 
within 1 kpc at late times, which drops significantly faster 
than the total sSFR of the galaxies. We interpret this as 
inside-out quenching during the post-blue nugget phase, fol¬ 
lowing the outside-in quenching during the pre-blue nugget 
phase. 

The central stellar mass density, a proxy for the bulge 
component, grows concurrently with the total stellar mass 
during the star-forming phase, typically at early redshifts 
and at low masses. At the high-mass end, there is a flat¬ 
tening of the Em*,i — relation, which is compatible with 
stellar mass growth in the outskirts for high-mass galaxies at 
late cosmic epochs. The simulations reproduce the observed 
— Mi, relation both at 2 ~ 2 and at 2 ~ 0. The mass 
growth in the outskirts at 1 — 5 7?e in the redshift range 
2 ~ 1 — 4 arises mostly from in-situ star formation and leads 
to a systematic gradual growth of stellar half-mass size: the 
simulated galaxies increase their sizes from sub-kpc scales 
to 2 — 3 kpc from 2 ~ 4 to 1, which is consistent with the 
sizes of early-type galaxies in the local universe. We empha¬ 
size that Re is evolving rather slowly with time compared 
to other quantities (such as the total stellar mass M*): Re 
grows by only a factor of ~ 4, while the stellar mass grows 
by a factor of ~ 100. 

Taking the comparison between observations and sim¬ 
ulations one step further, we have compared the profiles of 
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stellar mass, SFR, and sSFR at ^ = 2.2. Two important 
caveats in this comparison are that both the simulated and 
the observed samples are not mass- and/or SFR-complete 
(i.e., the observational sample selected mostly SFGs while 
the simulations did not apply any explicit selection on SFR) 
and that the SFR in the simulations and observations are ob¬ 
tained in different ways (counting newly born stars versus 
Ha luminosity). For the surface stellar mass density pro¬ 
files we find good agreement between the simulations and 
the observations of T15b, both for low-mass and high-mass 
galaxies, and both for the log-log profile shape and its zero- 
point normalization. For the SFR profile, the agreement in 
the profile shape is good, though the overall normalisation in 
the simulation is substantially lower by up to 1 dex, due to a 
certain overestimate of the SFR at very high redshifts in the 
simulations, as discussed in T15a. At z = 2.2, where most 
of the galaxies in our sample are in their post-blue nugget 
phase, we find for both observations and simulations a rising 
sSFR profile, i.e., the star-formation activity is reduced by 
an order of magnitude in the centre. This characteristic of 
inside-out quenching is stronger in the more massive galaxy 
sample (T15a). 

We conclude that the major compaction event and the 
subsequent quenching combine to a chain of events of ma¬ 
jor importance in the history of massive galaxies. In partic¬ 
ular, the high-redshift blue nuggets represent a key phase 
in galaxy evolution. The evolution of profiles of stellar and 
gas mass and of SFR tightly correlates with these charac¬ 
teristic phases of galaxy evolution, and can help us inter¬ 
pret them theoretically and identify them observationally. 
Our suggested phases can be directly observed. In the blue 
nugget phase, we predict a cuspy gas profile, whereas in the 
post-blue nugget phase, the gas profile is expected to have a 
hole, sometimes surrounded by a star-forming gas ring that 
quenches on a longer timescale. These gas signatures can be 
traced, for example, by IRAM Plateau de Bure and ALMA, 
which could resolve galaxies at 2 ~ 1 — 3 on sub-kpc scales. 
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Figure Al. Ratio of stellar mass to halo mass as a function of 
halo mass. The purple, green and blue point show our simulations 
at 2 = 2, 1.5 and 1.0, respectively. The large circles show the me¬ 
dian (at its 1 cr scatter) at these three redshifts. The horizontal 
black dashed line shows the cosmic baryon fraction of the Uni¬ 
verse. The solid lines show the (Behroozi et al. 2013a) relations. 
Our simulations lie a factor of 2 — 6 above the Behroozi et al. 
(2013a) relation at 2 ~ 1 — 2, but are in agreement with recent 
observations of SFGs at 0.5 ^ z < 2.6 (Burkert et al. 2015), which 
are shown as hexagons and their median as dashed line. A similar 
figure is shown in T15a. 
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APPENDIX A: STELLAR-TO-HALO MASS 
RELATION 

One of the key quantity for simulations to match is the stel¬ 
lar mass made within a given dark-matter halo. In this sec¬ 
tion, we compare the stellar-to-halo mass (M* — Mvir) rela¬ 
tion of the simulations with results from abundance match¬ 
ing (Conroy & Wechsler 2009; Moster et al. 2010, 2013; 
Behroozi et al. 2010, 2013a) and observed kinematics (Burk¬ 
ert et al. 2015). A similar analysis is presented in our com¬ 
panion paper (Tacchella et al. 2015a). 

Figure Al shows the ratio of stellar mass to halo mass as 
a function of halo mass for our simulations and the estimates 


from observational data (Behroozi et al. 2013a; Burkert et al. 
2015). At 2 = 2, we find for our simulated galaxies a median 
halo mass of log,Q Mvir = 11.5 ± 0.3 and a stellar to halo 
mass ratio of log,,, M*/Mvir = —1.4 ± 0.2. Via abundance 
matching, Behroozi et al. (2013a) found log,^Q M*/Mvir = 
—2.3 ± 0.2 at log,^Q Mvir = 11.5. However, Burkert et al. 
(2015) found logjg M*/Mvir = —1.6 ±0.3 at logjg Mvir = 
11.5, which is consistent with our simulations. 

Toward lower redshifts, the overall agreement between 
the simulations and observations improves. At 2 = 1, the 
simulated galaxies have a median halo mass of logj^g Mvir = 
11.7±0.3 and a stellar to halo mass ratio of logj^g M*/Mvir = 
—1.3 ± 0.2, which is again in agreement with the observa¬ 
tional estimates of Burkert et al. (2015), and slightly higher 
than the estimate of Behroozi et al. (2013a), who found 

logio M*/Mvir = -1.8 ± 0.2. 

We conclude that our simulations produce stellar to 
halo mass ratios that are in the ballpark of the values esti¬ 
mated from observations, and within the observational un¬ 
certainties, thus, indicating that our feedback model is ade¬ 
quate. 


APPENDIX B: PROFILES FOR ALL GALAXIES 

In Figures B1-B4 we plot the stellar mass, gas mass (total 
and cold gas), SFR, and sSFR profiles for all the 26 galaxies 
in our sample. The galaxies are order by decreasing total 
stellar mass at 2 = 2.0. 13 galaxies (V07, V09, Vll, V12, 
V20, V21, V22, V23, V26, V27, V29, V30, V34) show the 
clear signature of inside-out quenching, with a convergence 
of the stellar mass density in the inner regions. Four galaxies 
(V06, VOS, V13, V32) show only a weak signature of inside- 
out quenching, whereas nine galaxies (VOl, V02, V03, V15, 
VIO, V24, V25, V33, V34) show no such signature. 

This paper has been typeset from a TEX/IMEX file prepared by 
the author. 
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Figure Bl. Complete list of stellar mass (left), gas mass (middle-left), SFR (middle-right), and sSFR (right) profiles for all galaxies. 
We display the profiles at different redshifts bins [z = 6.0, 4.0, 3.0, 2.0,1.5,1.0) with different colours. The shaded regions correspond to 
the 1 a scatter of the profiles in a certain redshift bin. 
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Figure B2. Same as Figure Bl. 
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Figure B3. Same as Figure Bl. 
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Figure B4. Same as Figure Bl. 
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